Gradual eddy-wave crossover in superfluid turbulence. 
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We revise the theory of superfluid turbulence near the absolute zero of temperature and suggest 
a model with differential approximation for the energy fluxes in the fc-space, ennik) and £Kw(fc), 
carried, respectively, by the collective hydrodynamic (HD) motions of quantized vortex lines and 
by their individual uncorrelated motions known as Kelvin waves (KW). The model predicts energy 
spectra of the HD and the KW components of the system, EHnik) and Eawik), which experience 
a smooth crossover between different regimes of motion over a finite range of scales. For an exper- 
imentally relevant range of A = ln(^/a) {I is the mean intervortex separation and a is the vortex 
core radius) between 10 and 30 the total energy fiux e = eHD(fc) + eKw(fc) and the total energy 
spectrum £{k) — fHD(fc) + £'Kw(fc) are dominated by the HD motions for k < 2/£. In this re- 
gion £{k) follows the HD spectrum with constant energy flux e ~ Ehd =const.: £{k) oc k~^'^ for 
smaller k and tends to equipartition of the HD energy £{k) oc k^ for larger k. This bottleneck 
accumulation of the energy spectrum is milder than the one predicted before in [19] based on a 
model with sharp HD-KW transition. For A — 30, it results in prediction for the effective viscosity 
i/' ~ 0.002k (k is the circulation quantum) which is in a reasonable agreement with its experimental 
value in ''He low-temperature experiment ~ 0.003/s;. For k > 2/£, the energy spectrum is dominated 
by the KW component: almost flux-less KW component close to the thermodynamic equilibrium, 
£ « £kw «const at smaller k and the KW cascade spectrum £{k) -^ £Kw{k) ex: fc~^'^ at larger k. 



I. INTRODUCTION 

Liquid '^He and "^He at very low temperatures can be 
viewed as a perfect superfluid without any normal fluid 
present. Turbulence is a very common state for such su- 
perfluids, and it comprises one of the most interesting 
subjects in physics with exciting recent developments in 
turbulence; see, e.g. [1-16]. Superfluid turbulence con- 
sists of a tangle of quantized vortex lines [1-4, 17, 18]. 
How is superfluid turbulence related to usual hydrody- 
namical turbulence? On the one hand, at the scales 
greater than the mean distance between the inter-vortex 
separation distance one can expect the vortex discrete- 
ness to be unimportant and, therefore, superfluid and hy- 
drodynamical turbulence should have similar properties 
at these scales. This can be true, of course, if the su- 
perfluid vortex tangle is not completely random but po- 
larized and organized into vortex bundles which at large 
scales form similar motions as would continuous hydro- 
dynamic eddies. In turn, such a vortex polarization can 
be either introduced by the external forcing (yet to be un- 
derstood how), or it can occur due to a (yet to be found) 
self-organization mechanism. On the other hand, since 
there is no viscosity, the energy would cascade in super- 
fluid turbulence downscale without loss until it reaches to 
the small scales where the quantum discreteness of vor- 
ticity is important. It is believed that at this point the 
Kolmogorov-type (K41) eddy dominated cascade is re- 
placed by a cascade due to nonlincarly interacting Kelvin 
waves. Kelvin wave cascade takes energy further down- 
scale where it can be radiated away by phonons. 

Although the overall picture of superfluid turbulence 
described above seems quite reasonable, some important 
details of this picture are yet to be established. A par- 



ticularly interesting question is about the structure of 
the crossover between the eddy dominated and the wave 
dominated regions of the spectrum. As it was pointed 
out in our recent paper [19], the nonlinear transfer mech- 
anisms among weakly nonlinear Kelvin waves on discrete 
vortex lines is less efficient than the energy transfers due 
to the strongly nonlinear eddy-eddy interactions in con- 
tinuous fluids. This results is an energy cascade stag- 
nation at the crossover from the collective eddy domi- 
nated to the single-vortex wave dominated scales. The 
main message of paper [19] is that such a bottleneck phe- 
nomenon is robust and common for all the situations 
where the energy cascades experiences a continuous-to- 
discrete transition, and the details of particular mecha- 
nism of this transition are secondary. Indeed, most dis- 
crete physical processes are less efficient than their con- 
tinuous counterparts [26]. On the other hand, particular 
mechanisms of the continuous-to-discrete transition can 
obviously lead to different strengths of the bottleneck. 
The quantitative measure of the bottleneck is the value 
of rms vorticity, because this is the quantity which en- 
ters into the definition of the effective viscosity which 
is experimentally observable via measuring decay of the 
vortex line density [1, 7] (see below). 

Paper [19] considers the bottleneck mechanism under 
the simplest assumption that a sharp transition from the 
K41 eddy dominated cascade to the Kelvin wave weak 
turbulence occurs at the mean inter-vortex separation 
scale £. In this case, the vortex line reconnections pro- 
vide a mechanism to transfer the energy from the eddy to 
the wave motions, but their role for the energy cascade 
itself was neglected (even though it was noted that the 
bottleneck strength can be affected if this role was taken 
into account). On the other hand, paper [20] considered 



another extreme when the reconnections are the key pro- 
cess for the crossover cascade which was suggested to go 
through three different stages in a rather narrow range of 
scales of width A, - bundle-dominated, nearest neighbor 
and self reconnections. In spite of this rather unrealistic 
construction, the end result was still a bottleneck and 
reduction of the effective viscosity, though by a smaller 
factor than predicted in [19], A instead of A"''. 

In the present paper we neglect the role of reconnec- 
tions for the cascade process because, as we argued in 
[19], the reconnections are strongly inhibited within the 
polarized vortex bundles, and their occurrence is limited 
to the edges of these bundles. Since the volume in be- 
tween of the bundles is small compared to the volume in- 
side of these bundles, it seems natural to assume that the 
main contribution to the cascade will be due to nonlin- 
ear dynamics of non-rcconnccting vortex lines inside the 
vortex bundles, even though it is still possible that the 
reconnections can adjust the strength of the bottleneck, 
particularly if the K41 range is not too large and the tur- 
bulence polarization is reduced. The final answer about 
the role of reconnections should, of course, be sought in 
experimental and numerical data. In the present paper, 
we extend the analysis of [19] by taking into account the 
fact that Kelvin waves can be generated and play a role in 
the energy cascade at the scales greater than £, and that 
the transition from the eddy to the wave motions occurs 
over an extended range of scales rather than sharply. 



II. BOTTLENECK SCENARIO WITH SHARP 
CROSSOVERS 

A. Dimensional- and velocity-crossover scales 

When the "bottleneck" effect was first described 
in [19], it was assumed that the crossover from the eddy 
motions to the wave motions happens sharply at the scale 
~ £, i.e. mean intervortex distance, or in the fc-space at 
k ~ fcdim, where 



1 . 



(la) 



Subscript " dim " reminds that estimate (la) follows from 
the simplest possible dimensional reasoning. Besides, it 
roughly corresponds to an idea that at the scales larger 
than £ the vortex lines must be polarized and form bun- 
dles which would correspond, in a course grained sense, 
usual hydrodynamic eddies, while for k > kdim the vor- 
tex motions can be viewed as oscillations independently 
happening on individual vortex lines, i.e. as ID Kelvin 
waves. On the other hand, it was also remarked in [19] 
that the self-induced motion of the vortex line can get 
faster than its motion due to the collective interaction 
with the other vortices in the bundle already at the scales 
k <; fcvoi, where 



/Cvcl^ 



yVA. 



Subscript " vei " reminds that estimate (lb) follows 
from comparison of the self-induced velocity with cross- 
velocity induced at a given vortex line by nearby ^-distant 
vortex line. A more detailed discussion of this character- 
istic scale, which clarifies factor " 2 " under the square 
root, is given below in Appendix A 2. Thus, Kelvin waves 
can be expected to be present in some form already in 
the wave- vector range from k^oi to fcdim where they would 
coexist with the collective/eddy motions. It was argued, 
however, that due to their oscillatory character the waves 
would contribute much less into the cumulative motion 
of the vortex bundles in this scale range. 

Before going into details what occurs in the transition 
range 



fcvol r~^ f^ r~^ f^dh 



(Ic) 



it is worthwhile to reconsider some aspects of the problem 
under the simplest assumption that only one crossover 
scale is relevant. This is the subject of the current sec- 
tion. 



B. Bottleneck predictions in the presence of 
dimensional- and velocity-crossovers 

Here we reconsider the simplest scenarios of the eddy- 
wave transition with a sharp crossover with the only dif- 
ference from [19] that the crossover scale fc* is not neces- 
sarily at fcdim — 1/^ but lies somewhere in the range (Ic). 
To this end, we remind the Kozik-Svistunov [23] spec- 
trum of Kelvin waves, in the form suggested in [19]: 



£Kw(fc) ^ A{K'^e/£ 



8\1/5| 



-7/5 



(2) 



Here £Kw(fc) is the one-dimensional (in the fc-space) 
energy density of Kelvin waves, normalized such that 
-E'kw = J Snwik) dk is their total energy in unite vol- 
ume, K is the quantum circulation and e is the energy 
flux over scales. Parameters e and i arc mutually de- 
pendent, and their relation follows from the expression 
for the rms vorticity in the system of quantum fila- 
ments, -^/do^l^) ~ k£^'^. The later can be found for well 
developed turbulence when (|i-ij|^) is dominated by the 
classical-quantum crossover scale fc* . In the present case 



{\^\y 



>^ ^HDy"'/ ^f^ — "'^ ^K\V\"'*7 ' 



which, after substitution of (2), gives 
e ~ «VA5£i2^« 



(3) 



(4) 



Factor 2 in Eq. (3) follows from summation over vector 
indexes under assumption of isotropy of turbulent spec- 
tra. Equation (4) corresponds to the effective viscosity 



(lb) 



et 



K/A'^ikjf 



(5) 



Under the simplest assumption fc* ~ fcdim — 1/^ one gets 
the value, reported in [19]: v' ~ k/A^. For the sharp 
crossover at the velocity-crossover scale fcvoi — l/{iVA), 
where the self- induced velocity is of the order of the cross- 
induced velocity (see below), one gets v' ~ k/A. As we 
argued we expect that the true value of the crossover 
scale is somewhere in the region (Ic). Thus one expects 
that the true value of the bottleneck and corresponding 



is somewhere in between of the two extreme values 



kA-'' <v' <kA- 



(6) 



Note that formally the value i'' ~ kA^^ is the same as 
the one predicted by the Kozik and Svistunov approach 
based on reconnections [20]. 



C. Bottleneck at sharp amplitude-crossover 

Let us now check the consistency of the assumed in the 
previous section sharp crossover at some A* in the inter- 
val (Ic). To this goal let us evaluate the amplitude /i(fc,) 
of the Kelvin waves at this scale. Obviously, for consis- 
tency this amplitude must remain less than the intervor- 
tex separation h{k^) ^ £. This allows one to introduce 



the amplitude-crossover scale k^ 



at which 



h{k 



■amp J 



(7a) 



The estimate for h{k) can be obtained from the Hamil- 
tonian of Kelvin wave in the so-called local-induction ap- 
proximation (see, e.g., Eq. (5b) in our Rcf. [19]): 



h{k) ~ y'£Kw(fc)/A«;2fc3. 



(7b) 



Substituting here Snwik) from Eq. (2) and using Eqs. (4) 
and (7a) one gets estimate 



n^a 



#V, 



(7c) 



which is inside of the region (6). The subscript "amp" 
reminds that estimate (7c) follows from comparison of 
the wave amplitude with the intervortex distance. Factor 
" 2 " under the root is put by analogy with Eq. (lb) to 
ensure that fcvci < ^amp for any A. Assuming a sharp 
crossover at this scale one gets from Eq. (5): 



kA 



-11/3 



(8) 



which is, as expected, within the range (6). 

So, in order the waves amplitude to be less than the 
intervortex distance, the inequality k > fcamp rnust hold. 
The problem is to clarify what is going on in the interval 



k , < k < k 



(9) 



where formally computed [with the Kelvin-wave spec- 
trum (2)] wave amplitude h exceeds the intervortex dis- 
tance i^ which cannot physically happen. For example. 



at the scale fcvei one gets from Eq. (7b): h ~ A£ ^ i. 
On the other hand, the motions with k > fcvd cannot be 
considered as pure collective, because the cross-velocity 
(which is the influence of the motion of one vortex line in 
the place of another one) is smaller than the self-induced 
velocity (for more detailed discussion of this question, see 
Appendix A 2). 

Our scenario, is that in the interval (9) the growth of 
the wave amplitude on a particular vortex line would be 
arrested by the adjacent vortex lines in the bundle which 
would "get in the way". Speculations of similar type 
of Kozik and Svistunov [20] lead them to a suggestion 
that the hydrodynamic and the wave turbulence ranges 
are separated by the range of scales where the energy 
cascade is dominated by the vortex line reconnections. 
On the other hand, it was pointed out in [19] that, be- 
cause the vortex lines in turbulence must be polarized 
and organizes in bundles, the reconnection process must 
be suppressed and pushed to small volumes in between 
of the vortex bundles. Instead of a reconnection, one can 
expect a restriction of the wave motion of an individual 
vortex line when it grows in amplitude and tries to push 
close to the other vortex lines in the bundle. Naturally 
the growth of such a wave would get arrested at the am- 
plitude when the inter- vortex energy (which grows due to 
shortening of the distance to the considered vortex line) 
becomes equal to the vortex self-energy. 

This leads us to the following physical model of turbu- 
lence in the range (9). In the x-spacc, turbulence consists 
of vortex bundles with a fractal structure. Each vortex 
bundle which is made of denser sub-bundles, such that 
the mean separation of lines within the sub-bundle is 
^ i and the mean distance between the sub-bundles is 
3> £■ In turn, each sub-bundle consists of even denser 
sub-sub-bundles, etc. The density of vortex lines within 
a particular sub-bundle is such that at the scale of this 
sub-bundle the self-energy [which can be considered as 
the energy of the Kelvin waves £Kw(fc)] and inter- vortex 
energies [which should be associated with the hydrody- 
namic energy £iiu{k)] are balanced. This corresponds to 
condition in the k-space, 



^HDV^'J — ^KW\f^) : 



(10) 



in the range (9). For k > A:amp, Kelvin waves can prop- 
agate on an individual vortex line without approaching 
to (and being influenced by) the adjacent vortex lines. 
In the other words, the range k > fcamp is dominated 
by the Kelvin wave turbulence and the role of the eddy 
component will be clarified below. 

As we see, our corrected scenario which takes into ac- 
count that the eddy/wave crossover occurs over a finite 
range of scales predicts a bottleneck value which is in 
between of the values obtained by assuming sharp tran- 
sitions at the scales i and i^/A, respectively. 



III. 



FINITE CROSSOVER RANGE MODEL 



The goal of this section is to relax the simplified as- 
sumption that the bottleneek happens at some sharp 
crossover scale and to describe in a simple manner the 
transition regimes around the characteristic scales, in- 
troduced in the previous Section. The first step in this 
direction is to revise the differential approximation for 
the cascades of turbulent energy; this is done in the fol- 
lowing Subsection. 



A. Differential approximation for the cascades of 
turbulent energy 

The energy spectrum £{k, t) of isotropic turbulence can 
be described by the continuity equation 



d£{k,t) de{k,t) 



dt 



dk 



= 0, 



(11a) 



where e{k,t) is the turbulent energy flux over scales. In 
the stationary case this equation simplifies to the require- 
ment of the constancy of the energy flux in the so-called 
inertial interval, where both energy pumping and energy 
dissipation can be neglected: 



e{k) = e 



(lib) 



In order to describe a stationary spectrum £{k) one needs 
to know how e{k) depends on £{k). For simplicity in this 
paper we will use reasonably simple differential models, 
that describe the turbulent energy cascades of hydrody- 
namic (HD) and Kelvin wave (KW) turbulence at least 
qualitatively and sometimes even semi-quantitative. 

The first differential equation model for HD was first 
proposed by Leith in 1967 [21] and was recently studied 
in [22]: 



.»„W^4vinw«i^ 



(12) 



Here £HD(fc) is the energy flux carried by HD turbulence. 
For Eq. (13b), the factor | reproduces a numerical co- 
efficient that reasonably fits the experimentally observed 
value of the Kolmogorov constant. 

Generic HD spectrum with a constant energy flux was 
found in [22] as a solution to the equation e-aoik) = e = 
const: 



£-ar>\l^) — "^ 



24 £ /T\3/2- 

llfcll/2 + V^ 



2/3 



(13a) 



The large k range describes a thermalized part of the 
spectrum with equipartition of energy characterized by 
an effective temperature T, namely, T/2 of energy per 
a degree of freedom, thus, £k = Tk'^ j-Kp. At low fc, 
Eq. (13a) coincides with the K41 spectrum: 



£Ho(fc) = (24/11)'/^ £2/3 ^-5/3^ 



For Kelvin turbulence, the differential approximation 
model was suggested in [24]. In a way similar to Eq. (12), 
we suggest here a differential approximation for the en- 
ergy flux, carried by the Kelvin waves: 



m 



5 {kif£, 
7 



4 

KW 



(fc) d£^„{k) 



A5k7 



dk 



(14) 



Note that this form is slightly less general than the one 
of [24] because it does not take into account conservation 
of the waveaction. However, it is simpler which allows a 
more detailed analytical treatment. 

In the stationary case, equation £Kw(fc) = £ = const 
has the solution 



Ckw(,'^J 



(- 



1/5 



This solution changes from KW-spectrum for small k 



(15a) 



£^„{k) ~ A{K'^e/e 



8\ 1/5 ,^7/5 



(15b) 



to the thermodynamically equilibrium solution with 
equipartition of energy (Rayleigh-Jcans spectrum) 



£Kw{k) = T/irp, 



(15c) 



for large k. The factor — | in Eq. (14) is chosen such to 
reproduce in Eq. (15b) the numerical coefficient equal to 
unity. The actual value of this factor is still not estab- 
lished with a reasonable accuracy, see, e.g. [23]. 

When the eddy (HD) and KW turbulence coexist, both 
models should work together in such a manner that for 
small k the HD spectrum should be recovered, while for 
large k the KW spectrum should be the only one: 



£{k) = 



£^n{k), k <$:!/£, 

fKw(fc), k:^l/£ 



(16) 



A way to reach this physical requirement is presented in 
the following section. 



B. A unified model for the total eddy- wave energy 

flux 



A relatively simple model of turbulence with two types 
of motions, random eddies and Kelvin waves, is as fol- 
lows. The two types of motion coexist and interact in 
the extended crossover range in the following sense: 

a. The total turbulent energy density £{k) and the 
total energy fiux over scales, e{k), consist of two respec- 
tive parts: 



£{k) = £HD(fc)+£Kw(fc), 

e{k) = £HD(fc) +£Kw(fc) ; 



(17a) 
(17b) 



where energy fiuxes £hd(^) — £HD(fc)+£HD ^^'^ £Kw{k) — 
£Kw(fc) + £kw(^) have additional contributions e^-^ (k) 
and £kw(^) that originate from influence of KW on the 



(13b) HD-energy flux and vise versa. 



b. Continuity equations (11a) for the energy densities 
have to be supplemented by additional terms ±F{k), that 
describe energy exchange between two types of motion: 



d£ni}{k,t) denujk.t) 
dt dk 



dt 



dk 



-F{k,t), 
F{k,t), 



(18a) 
(18b) 



c. Cross-contributions to the energy fluxes, e-^^-^ (k) 
and EkwI^) ^^6 modeled in the linear approximation with 
respect of the influential energies (i.e. the HD energy 
influencing the KW flux and vice versa) : 

£„„ (fc) = PhdI^hd} d [£^„{k)/kl]/dk^ , (19a) 
e™(fc) - PkwI^kw} d [£„o(fc)/fc2]/dfc2 , (19b) 

with some wave-vector fc, which will be clarified later. 
Differential form of these contributions follows from phys- 
ical hypothesis that these terms should disappear (or be- 
came much smaller and can be neglected) when the in- 
fluential subsystem is in thermodynamical equilibrium, 
i.e. when £hd ot fc^ and £kw oc fc" =const. Functionals 
of the corresponding energies, !?...{...}, will be modeled 
by dimensional reasoning exactly in the way, how equa- 
tions (12) and (14) for the fluxes have been formulated. 
Resulting equations for A can be written in the form: 



-^hdI^hdJ 



L^HD y k CHD{k) , 

C/Kw(fct-) "'*^Kw("')'* 



(20a) 
(20b) 



where Chd is a dimcnsionlcss parameter and CKw{k£) is 
a dimensionless function of k£, that will be chosen below 
in Eq. (23). 

d. The energy distribution between the counterpart 
components depends only on k and for simplicity is as- 
sumed to be independent of the level of turbulence exci- 
tations: 

£nnik,t) = g{k£)£{k,t), (21a) 

£Kw(fc,i) = [l-g{k£)]£{k,t), (21b) 

were we introduced (only) fci?-depcndent blending func- 
tion g{k£) which will be explained below. 

e. Resulting model for the total energy flux e{k). 
Adding the two Eqs. (18) and using Eqs. (17), one 
yields the continuity Eq. (11a) in which £(k) is given by 
Eq. (17b). In such a way the unknown function F(k,t) 
disappears from the game. Together with Eqs. (12), (14), 
(17b), (19), (20) and (21) this flnally gives: 



e{k) = -{ 



y/k^^g{k£)£{k) 
5{k£fkl[l-g{k£)]^£{k)^ 



(22) 



7 



il^*^) 



A5k7 

gik£) 



1 - gjM) 

kl 



]}■ 



In the derivation of this equation we took 

Chb = -1/8 , Ckw(H) = -b{k£f/7K'' . (23) 

Only with this choice the resulting Eq. (22) for e(fc) is 
proportional to d[£{k)/k'^]/dk and e{k) vanishes in the 
thermodynamical equilibrium with £{k) ex fc^, as one 
should expect. 

Equation (23) contains yet unknown blending function 
g{k£) which will be discussed in the next section. 

C. Separation of the eddy and wave motions 

In order to find a qualitative form of the blending func- 
tion we consider a system of locally (in the vicinity of 
some point Tq) near-parallel vortex lines, separated by 
mean distance £ and supply them by index j. Notice 
that in principle the same vortex line can go far away and 
come close to Tq several times. To avoid this problem one 
should assign the same vortex line a different index j if 
it leaves (or enters) the ball of radius £\/X centered at 
rg. Each vortex line (with zero radius a) produces a ve- 
locity field Vj{r), which can be found by the Biot-Savart 
Law (Al). 

The total kinetic energy E ~ j ^^ {vi ■ Vj) can be 
divided into two parts, E = Ei + E2, where 



^1^JE(«')' S2^|E(^^-^^- 



«#j 



E 



V, ■ Vi 



(24) 

The same subdivision can be made also for the energy 
density in the (one-dimensional) fc-space, £{k) ~ £i{k) + 
£2(^)7 with two terms, that can be found via fc-Fourier 
components of the velocity fields Vj{k) in the way, similar 
to Eq. (24). Now our idea is as follows: energy £i{k) is 
defined by the form of the individual vortex lines, that 
is determined by the Kelvin waves, while energy £2{k) 
depends on correlations in the form of different vortices, 
that produce collective, hydrodynamic type of motions. 
Therefore £i{k) can be associated with the Kelvin wave 
energy, £i{k) =i> £^^{k), while £i{k) has to be associated 
with the hydrodynamic energy, £2{k) => £iiu{k). This 
allows one to conclude that 



g{k£)^[l+£i{k)/£2{k)y 



(25) 



The rest is technicalities presented in Appendix A, where 
we concluded that in practical calculations it is reason- 
able to use analytical form g{k£) of the blending function 



g{k£) = go [0.32 ln(A + 7.5) k£] , (26a) 



where 



go{k£) = 



{k£f exp(H) 
ATT{l + k£) 



(26b) 
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A 


/Cvcl ^ 

Eq. (lb) 


fccn-t 

Eq. (27) 


Eq. (7c) 


Numerics 


exlO"^ 
self-cons. 


10 


0.45 


2.3 


0.76 


24.4 


5.2 


30 


0.26 


1.8 


0.64 


23.6 


2.1 


10^ 


0.14 


1.4 


0.52 


22.8 


0.30 


10^ 


0.045 


0.96 


0.35 


22.5 


1.2x10"^ 



TABLE I: Comparison of the energy-, velocity-, amplitude- 
and flux-crossover scales for typical experimental values A — 
10 and 30 and unrealistically large values A = 100 and 10'^. 
Values of ka depend on e defined in a self-consistent way, 
explained in Sec. Ill D 3 



1. Comparison of various crossover scales 

With the proposed blending function we can introduce 
another cross-over scale ken^ by comparing HD and KW 
energies, at fcon^ they are equal: g{keni) = 1/2. As fol- 
lows from Eqs. (26), k^^i has very weak, logarithmical 
dependence on A: 



fccn^ ~ 6.64/ ln(A + 7.5), 



(27) 



presented in the second column in Table. I. 

First column of this table displays velocity-crossover 
scales, given by Eq. (lb): fcvoi^ — ■\/2/A, while in the col- 
umn 3 one finds the amplitude crossover fcamp^i Eq. (7c), 
scale at which the formally computed with the KW spec- 
trum amplitude of Kelvin waves reaches intervortex dis- 
tance. The fourth column of Table. I displays flux- 
crossover scale, kai at which the contribution of the eddy- 
and Kelvin-wave turbulence to the energy flux in the k- 
space become equal. This scale is introduced below in 
Sec. HID. 

Important message is that although theoretically, in 
the limit A — > oo (see Appendix), fcon <C fcamp, in the 
region 10 ;$ A < 30 these scales are close, moreover, due 
to numerical prefactors fcen > ^amp- Therefore the frac- 
tal structure of the vortex lines, described in Sec. II C, 
does not appear which allows us to use in the actual re- 
gion 10 ^ A ^ 30 the proposed simple blending function 
(26b). 



D. Turbulent energy spectra in the differential 
approximation 

1. Dimensionless representation 

At the end of the day we arc left with solving the ODE 
(22), where the blending function is g{k£), Eqs. (26b). At 
this stage of research it is reasonable to nondimensional- 
ize physical quantities introducing 



In particular, with this normalization the one- 
dimensional energy spectra HD and KW, Eqs. (13) and 
(15b), take the form 



eHo(x) = (24/ll)2/3e2/3^-5/3^ 



^HD V^) J- ) — ^ 



24 e 



lla;ii/2 



213/2 



2/3 



(29a) 
(29b) 
(29c) 



where T is non-dimensional temperature. And the 
ODE (22) to solve [with the boundary condition e(x) -^ 
eKw(a^) for a; ^ 1] becomes 



-{g\/2;"ff(a;)e(x) 

[l-g(a:)]V(x)}x 



(30) 



8^2 



X = kt , e{x) 



■£{x) 



(28) 



5 X x,^ 

— e(a;) 
ax 



Here we made the natural choice that the crossover scale 
/c* between two types of thermodynamic equilibrium is 
ken, the scale where energies of two types of motion are 
the same. 



2. A-dependence of the energy spectra 

An instructive solution e{x) with A = 30 and e — 
2.1 10~^ is shown in Fig. 1, left, as a (black) solid line. 
One sees that this solution for x < a;en — 1-83 follows the 
thermolized HD spectrum ei„:>{x,T) [given by Eq. (29b) 
with properly chosen T] shown as a dotted (cyan) line. 
An important observation is that the pseudo-thermolized 
part of the spectrum is very pronounced in the region 
X > 0.5 where it is very different from the K41 spec- 
tra of HD turbulence eHD(a;) oc x~^^^ shown as a dashed 
(blue) line. For x ^ xa the solution practically coin- 
cides with the pure KW spectrum eKw(2;), (29c) shown 
as a dashed (green) line. Important, that the crossover 
scale Xfi ~ 23.6, at which the total energy flux consists of 
50% of HD- and 50% of KW-fluxes, is much larger that 
Xcn — 1.83, at which a half of the total energy is carried 
by HD and half by KW motions. To make this evident 
we plotted in Fig. 1, right, the partial HD- and KW- 
energy fluxes vs. x. They become equal at a;fl, which for 
A = 100 and e = 3 10"'* is around 23.6. 

In the intermediate region x^n < x < xn the energy 
consists mostly of the KW energy, while the energy flux 
is carried mostly by the HD motions. Explanation to 
this observation is simple: as follows from Eq. (30) the 
HD motions are more effective (in factor ~ A^/x^^^) in 
support of the energy flux then the KW turbulence. Be- 
cause the main part of the energy flux is taken by the 
HD motions, the KW energy spectrum (and therefore 
the total one) is close to the flux-less KW-solution: ther- 
modynamic equilibrium (15c), £'kw =const. For x > Xfl, 
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FIG. 1: Left panel: Solid (black) line represents typical total dimensionless energy spectrum e{x) obtained by numerical 
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both the energy and the energy flux are carried by the 
KW motions. Therefore the total energy spectrum coin- 
cides with the KW cascade solution. For larger values of 
e the flux-crossover scale goes to the smaller values of k, 
see Fig. 2, remaining nevertheless larger than fccn- 

For smaller value of A qualitative behavior of £{x) re- 
mains the same, just different parts of the spectra (with 
larger values of self-consistent values of e) become less 
pronounced. 



3. Self- consistent estimate of the dimensionless energy flux 

Energy spectra, Figs. 1 and 2, which we obtained in the 
proposed differential approximation are quite similar to 



those suggested under the assumption of sharp crossover, 
see Fig. 1 in Ref. [19]: for small k they coincide with the 
HD spectrum, including the bottleneck part with (al- 
most) thermalized part £ oc k^ , while for large k the 
spectrum follows the KW spectrum oc fc~^/^. The only 
difference is that with sharp crossover the thermalized 
part of the HD-spectrum is matched with KW-spectrum 
at some k = kt,, while in the differential approximation 
with a smooth blending function there is essential inter- 
mediate region (about one decade) fccn < k < ka, with 
(almost) KW-thermalized part £ =const. This leads to 
an essential difference in the estimates of the vorticity 
(w^) and as a result in the estimates of the effective vis- 
cosity v' , a parameter that can be measured (implicitly) 
in experiments. Indeed, in the models with the sharp 



crossover one estimates (w^) in Eq. (3) as k^EHoik^,) and 
then equates £hd — i^kw, because in this model the HD- 
energy "transforms" into the KW-energy at the position 
of the sharp crossover. 

In the "continuous" model presented above one has to 
account for a wide region between ken and fcfl , where the 
flux is supported by the HD turbulence, while the en- 
ergy is dominated by the Kelvin waves. Therefore \oj'^) 
can be estimated similarly to Eq. (3) as 2 J k"^ £^^{k)dk ~ 
^on^HD(fcon), but now £HD(fccn) caunot be estimated based 
on the —7/5 Kozik-Svistunov spectrum because, as one 
sees in Fig. 2, £hd is much lower at this point. A s a re- 
sult, at the same energy flux the rms. vorticity a/(|wP) 
occurs to be essentially smaller then that in the sharp- 
crossover models and the effective viscosity is larger. In 
our approach the rms vorticity can be found more accu- 
rately by numerical calculation of the integral in Eq. (3) 
with spectrum £^i^{k) ~ g{k£)£{k), where the blending 
function is given by Eqs. (26): 



(M^)=2r 



Pg{k£)£{k)dk, 



(31a) 



where fcmin is the lower cutoff of the inertial interval. 
Using relation (|tt;p) = k^/£^ and normalization (28) one 
finds from Eq. (31a) in the limit k^in ^ 0: 



x^ g{x)e{x) dx 



(31b) 



Due to e dependence of the energy spectrum e, this rela- 
tion gives self-consistent a estimate of the dimensionless 
energy flux e, which is according to Eqs. (5) and (28) is 
nothing else, but v' / k. Resulting dependence e vs. A 
is shown in Fig. 2, right, as solid line. For convenience 
we approximate this dependence (in the actual interval 
A < 100) analytically: 



K 



.65 



103 + 45.8A + 1.98A2 



(32) 



shown in Fig. 2, right, as dashed line. Equation (32) 
reproduces numerical dependence e(A) with accuracy 
better than 1.5% for A < 50 and better than 8% for 
50 < A < 100. 

Notice, that found values of e = v' / n for A = 30 is 
0.021 which is quite close to the experimentally reported 
value v' ~ 0.003k in "^He experiments at low tempera- 
tures. A relationship between out model and experiments 
will be discussed below. 



E. 



Decay of quantum turbulence with the 
bottleneck energy accumulation 



Having in mind experiments with decaying supcrfluid 
turbulence, like the ones in [25], it is important to discuss 
how the bottleneck energy accumulation influences the 



decay of energy and vorticity in time. For this, we divide 
the total HD energy 



E„ 



into a sum of two parts: 



dk£^-a{k) 



Eno - -Bhd + -Eh 



(33a) 



(33b) 



the energy E^^^ associated with the K41 part of energy 
spectra £hd oc k~^'^, and the energy E^^^ associated with 
the thermodynamic equilibrium (TE) part of the spec- 
trum £hd fx k^. For our model: 



E^ 



TE 



dk£-ao{k) 



Cif£ C-HD [f^J 7 



(33c) 
(33d) 



fcTE 



where fcxE is the crossover scale between K41 and TE 
parts of the energy spectra corresponding to the posi- 
tion where £iio{k) is minimal. For A = 30, /cte ~ 1/^, 

K41 

see Fig. 1, left. The K41-energy, E^^,^ , is dominated by 
the outer region of the fc-space, k ^ kmm, while the TE- 
energy is determined by effectively the largest k ~ fcon 
of the HD motions (Fig. 1, left). In the experiment [25], 
^rnin is below 1 cm~^, which is less than fcen — 2/^ by one 
or two orders of magnitude. Then, the experiment [25] 

TE K41 

allows to estimate the ratio E^hd/^hd ^^ the proposed 
framework, and it varies from a few percents for small 
times to about 15-20 % at the latest times of the decay 
measurements. For us this means that with an accept- 

TE 

able accuracy one can neglect the contribution of E^^^ in 
Eq. (33b). 

Moreover, even when kept in Eq. (33b), the energy 

TE K41 

Sjjo does not appreciably affect on the decay rate of E^^ 
energy due to a large scale separation (fcmin <C k^^). The 

K41 

decay rate of E^^^ is determined by the energy flux e = 

K41 

—dE-^-^ I dt at the scale of the energy pumping, i.e. at the 
outer scale k = fcmin- The flux itself is proportional to 
£YiT,(kf^'^ ~ (i?jjj3 /krain)^^'^- For systcms with the time 
independent fcmin, as it is in [25], this gives the well known 
result for the late-time free-decaying HD turbulence: 



£^HD it) « t" 



(34a) 



and the time-evolution of the energy flux 

e{t) ex r^ . (34b) 

According to Eq. (28), £{t) = en^/i'^it) with the time- 
independent self-consistent dimensionless energy flux e, 
which depends only on A. This gives £{t) ex t^'^. There- 
fore, the vortex line density must decay in the standard 



1/f OC t-3/2 



(34c) 



in spite of the accumulation of energy E^^^ near the 
crossover scale fc„„. 

Notice, that energy E^j^{t) decays slower than 
_Ejjj5 (i) ex t^^. Indeed, in our model the dimensionless 
energy E^^f'/n'^ (cf. Eq. (28)) is dominated by the time 
independent scale Xen and, hence, by itself is time inde- 
pendent. Therefore, 



£;™(t)cx^-2g^r3/2 _ 



(34d) 



One concludes that E^^ energy is "decoupled" from 

K41 

the decay process of E^^ energy and does not affect the 
decay law (34c) of the vortex line density until to the very 
late stage of the decay, when the intervortex distance 
approaches the outer scale of turbulence and the entire 
model fails. 



Summary and Discussion 

In this paper, we revised the theory of the bottleneck 
crossover from the classical K41 cascade to the Kelvin 
wave cascade. In its previous form, transition from the 
eddy to the wave cascades was assumed to occur sharply 
at the scale i. The simple fact that the wave interactions 
are less efficient for the turbulent cascade than the hy- 
drodynamic eddies immediately yields prediction for the 
bottleneck accumulation near the crossover scale. How- 
ever, the bottleneck strength is rather sensitive to the 
details of the crossover region. In the present paper, we 
take into account that there exists a finite range where 
eddies and waves coexist and affect each other, making 
the crossover more gradual. As a result, the bottleneck 
in such a case is milder than in the model with the sharp 
crossover. To model the gradual transition range, we 
have employed a simplified turbulence model which is 
based on the differential approximation models of Leith 
type for the HD and KW components. Importantly, this 
model allows to make predictions for the realistic exper- 
imental values A in the range from 10 to 30, rather than 
making asymptotical predictions for the case A ^ oo. 
This appears to be important because, e.g., the asymp- 
totic theory gives fcamp ^ fcen whereas for A in the range 
from 10 to 30 we have fcamp ^ fcen- As a result, for the ex- 
perimentally important situations there is no range with 
equipartition of the eddy and the wave energies given by 
(10). For similar reasons, the theory of crossover [20] 
which fits three asymptotic ranges into a single decade of 
scales is rather unrealistic (leaving aside the issue about 
the role of reconncctions which we mentioned before). 

One may experience some problems trying to imagine 
any HD components at k(, > 27r, where the wavelength 
becomes larger than the intervortex distance, and even 
come to an idea that g{kt) must become zero sharply 
at fc£ — 2tt, or generally, k£ ^ 1. Our model is based 
on a (reasonable) hypothesis that the HD motions are 
identified with the coherent part of different vortex line 



motions. In such an approach there is no formal lim- 
itation for the value of k from above. At k£ > 1 the 
velocity produced by fc-distortion of a given vortex line 
in the position of another i?-separated line, which is the 
reason for correlations in their motions, decays exponen- 
tially with k£. That is why our blending function, which 
measures the fraction of the HD motions, decays expo- 
nentially with k£. The actual hypothesis in this place is 
that even when g{k£) <C 1 the nonlinear energy flux car- 
ried by the HD motions is governed by the same equations 
as for the pure HD motions when g{k£) ~ 1. We believe 
this is a step forward in comparison with the oversimpli- 
fied scenario of sharp crossover which, by the way, is a 
limiting case for our model where the blending function 
is just a unit-step function g{k) = 0(fcdim — k). 

The found value of e = i^'/k for A = 30 is 0.021 
which is quite close to the experimentally reported value 
I/' ~ 0.003k in "^He experiments at low temperatures. 
Having in mind that our model does not contain fit- 
ting parameters, optimists can consider this agreement 
as more than satisfactory. On the other hand, pessimists 
can consider any agreement with just one number as ac- 
cidental. Realists should recall that the suggested model 
is based on a hypothesis of blending function, which was 
estimated without taking into the account the vectorial 
structure of the velocity field and, moreover, includes 
very important (step-like) assumption about the pair- 
distribution function of the vortex positions which allows 
one to estimate sum (A7a) as integral (A7b). Our feel- 
ing is that these approximations do not affect the results 
too much and a good agreement between the model and 
experimental values of v' supports the suggested model. 
Notice that our model predicts not only the value of ly' 
but the entire energy spectrum, which consists of four 
parts: K41 HD energy spectrum with constant energy 
flux, £ ex fc"'"'/^, a HD equilibrium £ ex k^ , a KW equi- 
librium £ ~ const and a KW-spectrum with constant 
energy flux, £ ex k^'^'^. This very definite qualitative 
prediction calls for more detailed experimental and nu- 
merical study of the supcrfluid turbulence, which, as we 
believe, will support our model. 
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APPENDIX A: ESTIMATION OF THE HD-KW 
BLENDING FUNCTION 

1. Estimation of the velocity field, induced by the 
vortex distortion 

To estimate the HD-KW blending function gik£), given 
by Eq. (25) we consider a vortex line slightly distorted, 
say in y-direction, by a sinus with a small amplitude A, 
running along 2;-axis with the fc-vector k. Then M ~ 
(0, Akcos{kz + (f)), l)dz, where (f) is an arbitrary phase. 
Each hne produces velocity that can be found via Bio- 
Savart Law: 



approximately valid for any R: 



Vj{r) = 



47r 



+ 00 



d£. 



(Al) 



where Sj = r — Vj with r ^ being the radius- vector point- 
ing to the (.j - the length element along the j'-th vortex 
line and. The total mean density of the kinetic energy 
per unit mass is E = ^ (I^P); where V(r) = ^ ■ Vj{r) 
is the total velocity field, and one understands (. . .) as 
the averaging with respect of the (random) vortex-line 
positions. 

Using Eq. (Al) we can find the resulting velocity (for 
simplicity) in the ZK-plane at distance R from the line. 
This velocity has only one component, Vx, and it is given 
by: 

K f°° R-Asa\{kz + (i)) + Akzcos(kz + (p) , , , „, 
^^ = ^ — ^ —T^dz.(k2) 



iz'^ + [R-Aiiin{kz + (t))f\ 



3/2 



Computing the contribution to the amplitude of the ve- 
locity variations proportional to the distortion amplitude 
A, one finds magnitude [i.e. factor in front of cos(fcz + (/>)] 
of the velocity fluctuations: 



5v{R) 



Sv{R) 



kA 



■kRx 



2ttR^ 
ikR[K2{kR)-Ko{kR)]-Ki{kR)V 

kA 



2.i?^V^+2'^^^P(-'^) 



(A3a) 



(A3b) 



where K„ are modified Bessel functions of the second 
kind. Interesting, that a simple interpolation (A3b) re- 
produces the exact result (A3a) with an accuracy better 
than 2%. 

Equations (A3) assume that the distance R to the vor- 
tex line is large enough to be able to neglect finiteness of 
the core radius a. In the opposite limit, i? <C a, one can 
use an equation for the self-induced velocity, in which the 
core radius is hidden in the parameter A — hi{£/a): 



(5ws 



AAk'^K/Air . 



(A3c) 



Considering this equation formally as a limit for the ve- 
locity when i? -^ 0, we suggest an interpolation formula 



where subscript " k " reminds that this velocity is induced 
by vortex line, distorted by sin(fcr). 



2. Velocity crossover scale 

A way to estimate the eddy-wave crossover scale is to 
compare the self-induced velocity filed with the veloc- 
ity field produced by a different neighboring vortex line. 
For this goal we consider simple geometry with two sin- 
disturbed parallel vortex lines separated by distance i 
They "start to feel" each other when 6v ~ 6v , i.e.: 



Sv ~ Sv 



<=;> 



v/4 + 27r fcvci^ 



f^vfil i- 



A(fc™i^)' . (A4) 



An approximate solution to this equation (with an accu- 
racy better than 10% for A > 10) is given by Eq. (lb). 
One can think that for the scales less than kvd, there 
is the cumulative effect and hydrodynamic-like behavior. 
For the scales larger than fcvoi the vortices does not seem 
to feel each other much, and the separate vortex line be- 
havior is important, hence the Kelvin waves. Corrections 
to this simplest viewpoint will be discussed below. 



3. Estimate of the blending function 

Having at hand an estimate of the velocity filed (A3d) 
we can use Eq. (25) to find the blending function g{ki). 

a. The first step is to obtain the kinetic energy 
density Eii(k) (per unit length), produced by one sin- 
distorted (with wave-vector k) vortex line: 



Eiiik) = 



Po 



Vk{R,z)\ ) dxdy (A5a) 



(AKAfc)2 p(l+^p/2)e-2p ^ .^^^. 

Po- — s ,^ A 9X9 Pdp (A5b) 



Stt 7o (2 + Ap2)2 

A 2 

PQ-—(AKk) , for A> 1 

Stt ^ ' 



(A5c) 



Here the fluid density po should be distinguished from 
the dimensionless radius p = kR, R = \/ x^ -V tp . In 
Eq. (A5a) we should substitute Vk{R, z) ~ Vk{R) cos{kz+ 
(j)) from Eq. (A3d) integrate over (x, j/)-plane orthogonal 
to the mean vortex line directed along 2;, and average 
(cos^ {kz + 0)) = i along z. Considering the integral in 
dimensionless polar coordinates p and ip after the free 
integration over i^, we get Eq. (A5b), in which integral 
can be analytically taken in the limit of large A with the 
result (A5c). This formula can also be obtained directly 
from the Kelvin- wave Hamiltonian in the so-called local- 
induced approximation, see e.g. Eq. (5b) in [19], where 
one expands -^/l + \dw{z) / dz\'^ over the line distortion 
w(z) = ^sin(fcz). 
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b. The second step is to find cross-kinetic energy 
density Eij(k, £ij) (per unite length), proportional to the 
product of the velocities produced by two sin-distorted 
vortex lines, separated by iij. 



OO /'OO 



Eij{k,£ij) ^ po dxj dy{\vk{Ri,z)vk{R2,z)\) 

J — OO J — OO 

OO ^ij/2 

= y (^)' j^- j SynPi)F{p2) 6-"+, (A6a) 

(A6b) 



yi + 7TP/2 

F\P) = o , A 9 ' P+= Pl+ P2 



(.ij == ke.ij, Pi = A:ya;2 + j/2 



2 + Ap2 ' 

Here x = kx, y = ky, 
and p2 = k^ x^ + (j/ ~ ^ij)^- The change in the limits 
of integration is due to the (average) symmetry of the 
integrand relative to the plane y — £ij/1. The main 
contribution to the integral comes from the region with 
X ^ iij, which allows to simplify the integrand replac- 
ing p2 => tij — y + x^/{2lij). This clarifies the leading 
dependence Eij{k,£ij) on £ij as 



Eij{k,£ij) ex exp(— fc^ij) . 



(A6c) 



Algebraic improvement of this estimate cane be ob- 
tained by neglecting the (x, y)-dependence of F{p2) in 
Eq. (A6a). This procedure can be justified asymp- 
totically for k£ij 3> 1. The resulting dependence 
Eij{k,£ij) ex F{k£ij)cxp{—k £ij) approximates the in- 
tegral Eq. (A6a) with an accuracy of 10% for k£ij > 2. 
For smaller k£ij one has to account for £ij dependence 

of the remained integral /_^ . . .dy ~ (a/ A) ln(fc^ij -I- b) 
with a and b weakly dependent on A. As a result we 
approximate (with accuracy about 25%) the k£ij depen- 
dence of the energy Eij as follows: 



Eij{k,£i,) 



A 2 

po-^ (Auk) aln{k£ij+b)x 



8tt 

Vl + 7rfc4j-/2 
2 + A(Hij)2 



exp{~k£ij) . (A6d) 



In typical experiments with supcrfluid "^He and ^Hc the 
value of A varies from 10 to 30. Therefore we present here 



parameters a « 1.37, b ss 1.25 for A = 10 and a sa 2.37, 
&« 1.17 for A = 30. 

c. The third step is to find the relative total cross- 
energy of all j 7^ 1 pairs 



Rik£) = J2 



j¥1 
47r 
(H)2 



Ei,{k,£ij) 
En{k) 



(A7a) 



''"'^'^Wr- '""' 



which is estimated in Eq. (A7b) in continuous approx- 
imation assuming that j-lines are randomly distributed 
around line i = 1 with mean density l/£'^ for £ij ex- 
ceeding, value £. Generally speaking, the integral (A7b) 
should contain probability function of vortex separations 
V{£ij). Having no reasonable model for it, we choose a 
simple step function V{x) — for x > 1. 

Blending function g{k£), defined by Eq. (25), is related 
to the ratio (A7) as follows: 



g{k£) ^ R{k£) / [1 + R{k£)] , 



(A8) 



Taking oversimplified representation (A6c) of inte- 
gral (A6a) one computes go{k£), see Eq. (26b): 



9oik£) 



{k£f exp(H) 



47r(l + k£) 



(A9) 



Comparing A-independent function go{k£) with gnumik£), 
which depends also on A, we have improved its analytical 
representation by introducing A-dependent rescaling of 
the argument, see Eq. (26a): 



g{k£) = go[0.32\n{A + 7.5) k£] 



(AlO) 



The resulting function g{k£) gives a very reasonable ap- 
proximation to the results of "exact" numerical calcu- 
lation of gnum{k£). Therefore in practical calculations 
we will use analytical form (26) of the blending function 
gik£). 
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